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ABSTRACT 

Context. Our understanding of the formation and early evolution of globular clusters (GCs) has been totally overthrown with the 
discovery of the peculiar chemical properties of their long-lived host stars. 

Aims. As a consequence, the interpretation of the observed color-magnitude diagrams and of the properties of the GC stellar popula¬ 
tions requires the use of stellar models computed with relevant chemical compositions. 

Methods. We present a grid of 224 stellar evolution models for low-mass stars with initial masses between 0.3 and 1.0 M© and initial 
helium mass fraction between 0.248 and 0.8 computed for [Fe/H]=-1.75 with the stellar evolution code STAREVOF. This grid is 
made available to the community. 

Results. We explore the implications of the assumed initial chemical distribution for the main properties of the stellar models: 
evolution paths in the Hertzsprung-Russel diagram (HRD), duration and characteristics of the main evolutionary phases, and the 
chemical nature of the white dwarf remnants. We also provide the ranges in initial stellar mass and helium content of the stars that 
populate the different regions of the HRD at the ages of 10 and 13.4 Gyr, which are typical for Galactic GCs. 
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1. Multiple stellar populations in globular clusters: 
Definitions and possible origins 

Globular clusters (GCs) are described in classical textbooks as 
old, massive, tightly bound systems hosting single stellar popu¬ 
lations, i.e., coeval stars sharing the same initial chemical com¬ 
position. These objects are extensively used as ideal benchmarks 
for stellar evolution theory, and their estimated absolute and rel¬ 
ative ages provide markers for the formation and evolution of 
their host galaxies and their substructures fe.g. jKruiissenll2014b . 

This idealistic picture has been shaken with the discov¬ 
ery of multiple stellar populations in individual Galactic GCs 
thanks to high-resolution observations either with the Hubble 
Space Telescope or from the ground. Deep photometric stud¬ 
ies have revealed the presence of multiple sequences or spreads 
at the turn-off and/or along the subgiant branch Jn the color- 
magnitude diagrams of several GCs (e . g.. iBedin et al.l 


D’Antona et alTl2Q05t iPiotto et al J 12007b Villanova et al.1 


Mil one et all 2008 ~l2010l: Piottof 20081120091: iDalessandro et al] 

201 ll: iKing et al.ll2012l: iMilone et al.ll2013l) . In addition, spec 


20041: 

2007: 


troscopic analysis is now commonly used to differentiate be¬ 
tween so-called first- and second-generation GC stars (here¬ 
after 1G and 2G) based on their positions _along the^ well- 
documented O-Na anti-correlation (e.gJPrantzos & Charbonnell 
120061: ICarretta et al.l l2009bl: ICarrettal 120131) . This remarkable 
chemical pattern, which has been observed so far only among 
GC stars (see, e.g.. iBragaglia et al.1l20l4 iMacLean et al.ll2015[ 


and references thereinfl is advocated to be the main feature that 
distinguishes GCs from other stellar systems (e.g.. ICarretta et al.l 
2010 ). 

Together with the C-N and Mg-Al anti-correlations, it is in¬ 
terpreted as the signature of hydrogen-burning at high temper¬ 
ature in massive, fast-evolving 1G GC stars (hereafter referred 
to as polluters) whose ejecta (depleted in carbon, oxygen, and 
magnesium, and enriched in nitrogen, sodium, aluminum, and 
helium) mixed to various degrees with original intra -clus ter ma¬ 
terial and gave birth to 2G stars (e.g JPrantzos et al.1l2007b . 

Different scenarios of secondary star formation in¬ 
voke different types of 1G polluters, namely fast rotat¬ 
ing massive stars (FRMS) w ith initial masses above 25 
Me (iMaeder & Mevnet 120061: iPrantzos & Charbonnell 120061: 


iDecressin et alJ 120073114 [Krause et all 12017lt . massive asymp¬ 
totic giant bra nch (AGB) stars with init i al masses between ~ 
6 and 11 AC, (lYentura et al.lI2001L EoTl iD’Ercole et al.lI2010L 
1201 ILI20I2I: I Ventura & D’Antonal201~ ])7and supermassive stars 
with initial masses around 10 4 Me (iDenissenkov & Hartwickl 
120141) . In some cases the possible contribution of massive bi¬ 
nary stars (Ide Mink et al. 20091: Tzzard et al.ll2013l) . of FRMS 
paired with AGB stars dsills & Glebbeek l2010b. and of FRMS 
paired with high-mass interactive binaries feastian et al.ll2013bl: 
ICassisi & Salarisll2014l) are also considered. 

Several studies have suggested that the concept of stellar 
generations may not be appropriate to describe the long-lived 
stellar populations we are observing today in GCs, as this implies 
discontinuous star-forming events or long-term star formation. 


E-mail: william.chantereau @ unige.eh 


1 This feature has also recently been discovered among a small frac¬ 

tion of halo field s t ars (up to 2.5 % of the star sample); see, e. g., 

ICarretta et alJ (l2010l) : lMartell & Grebe! (l2010l) : lRamfrez et al.l (120121) . 
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In particular the need for multiple star-forming episodes occur¬ 
ring 30 - 100 Myr after the formation of the 1G, as required by 
the AGB scenario, is in conflict with the star formation history 
derived for young massi ve clusters in the Milky Way and in Lo¬ 
cal G roup galaxies (e.g jBastian et al .112013 at ICabrera-Ziri et al.l 
120141) . An additional difficulty of this delayed star formation 
scenario lies in the need to re-accrete Galactic gas with exactly 
the same che mical composition as that of the original proto-GC 
material (e.g.. ID’Ercole et al1l201ll I20I2I) after several tens of 
Myrs and in any case after the expulsion of the original gas and 
of the SNejdecta by SNe explosions or dark-remnant activation 
dKrause et al.ll2013l) . The probability that this has happened for 
all individual GCs during their Galactic journey is actually very 
small. 


At the opposite end in the FRMS scenario, the forma¬ 
tion of the sodium-rich, oxygen-poor low-mass stars (so-called 
2G stars) is predicted to be completed before fast gas expul- 
sion within ~ 3.5 - 8 Myr after the formation of the polluters 
dKrause et al .1120131) . This timescale is compatible with the fact 
that very massive star clusters with ages of o nly a few My rs in 
nearby galaxies appear to be devoid of gas dBastian & S traded 
120141: iBastian et al.ll2014l) . It is also shorter than the pre-main 


sequence evolution of the long-lived low-mass stars we observe 
today in GCfl; low-mass stars formed out of original gas and 
those formed out of polluted material in the immediate vicin¬ 
ity of the FRMS should therefore have arrived on the zero age 
main sequence with very small time delay compared to the typ¬ 
ical age of a GC of ~ 13.4 Gyr. It was also proposed that 
the ejecta of FRMS and/or high-mass interacting binaries were 
actually accreted by low -mass pre-main sequence stars of the 
same generation (jBastian et~aLll2013bl: see also Scenario II of 
iPrantzos & Charbonnell2006lfl in which case the concept of dis¬ 
tinct stellar generations is also misleading. Last but not least, 
in the FRMS polluter model a large percentage of (if not all) 
the stars exhibiting low sodium abundancejand that are there¬ 
fore countecias 1G members; see iPrantzos & Charbonnelll20061: 
ICarretta et al.ll2QQ9bh could have formed from pristine material 
mixed with th eH -burning ashes of the fast-evolving 1G polluters 
dCharbonnel et al.1120141) : in this case GCs would be formally de¬ 
void of 1G low-mass stars, their hosts being actually secondary 
star formation products with chemical properties spanning the 
whole range of the observed O-Na anti-correlations. 


Consequently, to avoid conceptual as well as semantic ambi¬ 
guities, we adopt the denominations “first population " (hereafter 
IP) for the stars that were born with the original composition of 
the proto-GC, and “secondpopulation " ( 2P ) for the stars that are 
born from a mixture between original material and the ejecta of 
short-lived IP polluter stars. 


2 The pre-main sequence lifetimes of our 0.8 M© models computed at 
the metallicity of NGC 6752 with initial helium mass fractions of 0.248, 
0.4, and 0.8 are ~ 1.045xl0 8 ,4.32xl0 7 , and2.87xl0 2 3 * * 6 yrs, respectively 
(see Sect. 0. 

3 A possible difficulty is that in this case a very short accretion phase (~ 

1-2 Myr) is necessary to acco unt for helium enrichment and for the ob¬ 

served Li-Na anti-correlation dD’Antona et al.|[2014l ; ICassisi & Sa laris 

120141 : ISalaris & Cassisill2014l) . 


2. The importance of quantifying helium variations 
among GC stars 

2.1. Dependence of the initial helium spread on the polluter 
scenario 


Despite the many differences between the scenarios that have 
been developed to explain the abundance properties of GCs, it is 
widely accepted that the Na-enriched ( 2P ) GC stars have started 
their life with a higher helium content than their IP counter¬ 
parts (this is actually obvious as helium is the main product of 
hydrogen burning). This has been confirmed by direct spectro¬ 
scopic measurements of non-local thermodynamic equilibrium 
He abundances for a sub set o f blue horizontal branch stars in one 
GC (NGC 2808: lMa.rino eta.U20l4 . However the initial helium 
distribution among 2P stars is not well known yet, as predictions 
for the extent of helium enrichment along the sodium distribu¬ 
tion in 2P stars strongly depend on the nature of the invoked IP 
polluters. 

In the case of the AGB scenario, the quantity of helium 
released by the polluters is set during the second dredge-up 
episode on the early AGB; it amounts to a maximum of ~ 0.36 
- 0.38 in mass fractio n, independe ntly of the initial mass of the 
AGB progenitor (e.g jDohertv et al.ll2014l) . On the other hand, 
the sodium content of the ejecta of the polluter results from 
the competition between hot-bottom-burning and third dredge- 
up episodes that may occur during the later thermal-pulse AGB 
(TP-AGB) phase; therefore, it strongly dep ends o n the initial 
mass of the polluter dForestini & Charbonnelll 19971: ISiessll2007l: 
ID’Ercole et al.l l2010t lYentura et al.l 120131) . Consequently, 2P 
stars spanning a large range of Na abundances are all expected to 
be born with very similar helium contents (maximum of 0.36 - 
0.38 in mass fraction if no dilution with the ISM matter is taken 
into account compared with ~ 0.248 for IP stars; see § 13.2.11) . 


In contrast, in the FRMS scenario the enrichment in sodium 
and helium of the polluters’ ejecta results from central hydrogen¬ 
burning on the main sequence. Therefore this model predicts 
broad and correlated spreads in both sodium and helium among 
2P stars at birth. In the case of NGC 6752, for example, 2P 
stars are expected to be born with initial heliu m mass fractions 
ranging between 0.248 (as for the IP) and 0.8 dDecressin et al.l 
l2007bh . Such a large helium dispersion correlated with the 
sodium spread provides a straightforwar d explanation to the lack 
of sodium-rich 2 P AGB stars in this GC dCharbonnel et al.l2013l: 
observations bv ICampbell et al.ll2013l) due to the impact of the 
initial helium content on the stellar lifetime and evolution path 
in the Hertzsprung-Russel diagram. Simultaneously, the FRMS 
scenario predicts that today (i.e., at a typical GC age of 13.4 
Gyr), 95% of the low-mass stars lying two magnitudes below 
the turn-off should be born with a helium mass fraction Yi n i be¬ 
tween 0.248 and 0.4, an d only 5 % with Yi n i higher than 0.4 
dCharbonnel et al.l 120131) . Additionally (see Sect. [5]), no star 
born with Yi n i higher than 0.4 is expected to lie today on the 
horizontal branch, in agreement with current interpretations of 
GC horizontal branch (HB) morphologies in conn ecti on with th e 
so-called s econd parame ter probl em (e .g., D’Antona et al. 2002; 

D’ Antona & Caloi|l2004tlCaloi & D’ Antonall2007l:lGratton et al] 

20101:iDotter et aklboiOh. 


2.2. Objectives of the present series of papers 

The aim of this series of papers is to investigate the implications 
of the initial helium spread predicted by the FRMS scenario in 
different parts of the GC color-magnitude diagram, and at vari- 
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ous ages along GC evolution. Predictions will be compared to 
those of the AGB scenario and to available relevant observations, 
and future observational tests will be suggested. 

Here (Paper I) we present a grid of low-mass, low-metallicity 
stellar models that will be extensively used for different pur¬ 
poses in the series. These models are computed with a [Fe/H] 
value of -1.75 close to that of the well-studied GC NGC 6752 
and for a wide range of initial helium abundance, taking into ac¬ 
count the corresponding variations expected within the FRMS 
model for the elements involved in the CNO cycle and the NeNa 
and MgAl chains (see § [3] for a detailed description of the ini¬ 
tial chemical compositions and the input physics of the stellar 
models). We discuss the impact of the initial chemical compo¬ 
sition on the characteristics, lifetimes, evolution behavior, and 
fate of GC low-mass stars; we first focus on the case of a 0.8 M© 
star (§ 0 ), and then present the dependence on the initial stellar 
mass (§0. Although the influence of initial helium abundance 


on the stellar properties has already been extensively described 

in the literature (e.g., Iben & Faulkn^ 1968s; 1 Iben & Rood 

1961 

Demarque et al. 1975; Sweigart 1978; D’Antona et al. 

2005; 

Caloi & D’Antona 20071: Salaris & Cassis 

2005) 

Maeder 

2009; 

Pietrinferni et al. 2009); Sbordone et al.| 

2011 ) 

Valcarce et al. 

2012; Cassisi et al. 201j|; Kippenhahn et a] 

. 2013|), to our knowl- 


edge this is the first investigation of the impact of helium content 
higher than ~ 0.42 on the whole evolution (from the pre-main se¬ 
quence to the end ojf the AGB phase) at low [Fe/H], and for such 
a broad mass rang 


3. Grid of IP and 2P low-mass stellar models at 
[Fe/H] = -1.75 

3.1. Generalities 

We present a grid of standard stellar models (i.e., with no atomic 
diffusion or rotation, and no overshooting) with initial masses 
between 0.3 and 1.0 M© (mass steps of 0.05 M 0 ) for [Fe/H] close 
to that of NGC 6752, and with Yi n i ranging between 0.248 and 
0.8 and the abundances of the elements included in the CNO, 
NeNa, and MgAl chains that change accordingly (see § 13.21 for 
details on the initial chemical mixture assumed in the computa¬ 
tions). The models are computed from the pre-main sequence 
up to the tip of the TP-AGB or the pl anetary neb ula phase 
with the evolution code STAREVOL (e. g.,[Palacios etliill20061) . 
Mass loss is treated assuming the iReimersI ( 19751) formula (77 
parameter of 0.5) on the red giant branch (RGB) and d uring 
central helium-burning, and the IVassiliadis & Wood (1199 3) pre¬ 
scription on the AGB. The basic input physics (nuclear reac¬ 
tion rates, equation of state, opacities, mass-loss prescriptions) 
is the same as in lLagarde et al. ( 2012b . except for the mixing 
lengt h parameter (1.75 here, as assumed in the FRMS models of 
IDecressin et al]l2007bl instead of 1.6). 

3.2. Initial chemical composition 
3.2.1. IP models 

The initial composition of IP stellar models must reflect that 
of the original intra-cluster gas. We assume an initial helium 
mass fraction Yi n i of 0.248 (i.e., slightly higher than the pri¬ 
mordial he lium a bundance 0.2463+0.0003 predicted by standard 
BBN bv ICoc et~aPl2013l using the new value of f \^ determined 

4 iD’Antona et al.l ( 20101) have investigated the effect of a Y up to 0.8, 
but only on the horizontal branch morphology. 


bv iPlanck Collaboration et afll2013l) . In order to be consistent 
with IDecressin et a l. ( 2007b) we adopt a scaled-solar abundance 
composition dGrevesse & Noelsl Il993h . a heavy element mass 
fraction Z of 5 xlO -4 , and an alpha-enhancement of [cr/Fe] = 
+ 0.3. This corresponds to [Fe/H] = -1.75, close to the value of 
-1.56 derived for NG C 6752 (e.g.. ICarretta et al.lI2007L l2Q09at 
IVillanova et al.ll200~9l) . 

3.2.2. 2P models 

Following iPrantzos & Charbonnell ( 2006 ) and IDecressin et al.l 
(2007a, b), we assume that the initial chemical composition of 2P 
stars consists of a time-dependent mixture between a) the intra¬ 
cluster material of the original composition and b) the mate¬ 
rial processed to various degrees through central H-burning and 
progressively ejected by the FRMS during their mechanically- 
driven slow wind phase (i.e., along the main sequence and the 
Dr-limit phase). As FRMS release the first material of origi¬ 
nal GC composition at the very beginning of their evolution, and 
later material richer in H-burning ashes, one expects 2P stars 
with Na and He contents similar to that of IP stars to form first, 
and then 2P stars with various deg rees of Na and He enrichment 
(see, e.g JCharbonnel et afll2014l) . 

To account for this general behavior, we assume that the ini¬ 
tial mass fraction of a given isotope i in 2P stars varies from star 
to star according to 

X 2P = x wind (1 _ a( ) + X' CM a t , (1) 

with the various terms described below. 

- The original intra-cluster medium value X ICM is equal to Xj F 
assumed for IP stars (§ 13.2.11) . 

- For the mass fraction of the isotope i in the massive star 
ejecta, X wmd , we use the time-dependent predictions along 
the main sequence and the Dr-limit phase of the 60 M© 60rC 
model FRMS model of IDecressin et al.l ( 2007b ) computed 
with the Set C of nuclear reactions (see their Table 2, and the 
discussion in § 13.2.31) . In order to cover the highest helium 
enrichment, we also take into account the part of the ejecta 
of their 120 M© 120rC FRMS model that correspond to a he¬ 
lium mass fraction higher than 0.88 (i.e., the maximum value 
for the 60 M© star; see Fig.Q]). 

- The dilution factor a t varies during the lifetime of individual 
IP polluters between 0.99 (ejecta strongly diluted with the 
surrounding intra-cluster gas) and a minimum value ami n of 
0.2 that is constrained by the Na-Li anti-correl ation ob served 
in NG C 6752 turn-off stars (observati ons bv iPasquini et al.l 
12005b see, e.g.. IDecressin et al.1l2007bl) . 

The resulting range for the initial abundances of C, N, O, 
Na, Mg, and Al of our 2P stellar models is shown in Fig. [I] as a 
function of the initial mass fraction of helium. For each stellar 
mass between 0.3 and 1.0 M© (mass steps of -0.05 M©, i.e., 
14 initial masses), we compute 16 models with initial helium 
abundance varying between 0.248 and 0.8 (steps < 0.05, actual 
Yi n i values indicated by the black points in Fig. []]). Therefore, 
our grid consists of 224 stellar models in total. 

3.2.3. Comparison with spectroscopic abundances in 
NGC 6752 

Within the FRMS scenario, 2P low-mass stars are expected to 
form in the immediate vicinity of the massive star polluters, 
which might actually all evolve in different environments (e.g., 
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Fig. 1 . Theoretical distribution for the initial carbon, nitrogen, oxygen, sodium, magnesium, and aluminum abundances ( [X/Fe] = log(X/Fe)*- 
log(X/Fe) 0 ) for [C/N] and for 12 C/ 13 C (initial value of 90) as a function of the initial helium mass fraction adopted for our 2P stars models (black 
line; dots indicate the 16 values adopted for the model computations of each stellar mass between 0.3 and 1 M 0 ). This accounts for dilution 
between original IP material and the ejecta of the 60 and 120 M© FRMS models (red and blue lines, respectively; iDecressin et alJl20 07bh as 
described in Eq. 1 ( see text). The mage nta arrows represent the va riations observed in NGC 6752 for [C/Fe] and fN/Fel (ICarretta et al.ll2005h , 
[O/Fe] and [Na/Fe] (ICarretta et al .1120071) . and [Mg/Fe] and [Al/Fe] (ICarretta et al.ll2012h 

in terms of pristine gas available), depending on their orbit dKrause et al.l2013l) . Individual FRMS polluters of similar initial 
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Fig. 2. (Top) Evolution in the HRD of the 0.8 M© models for 16 
different Y ini between 0.248 and 0.8 (Z = 5 x 10 -4 ). Colors depict the 
evolutionary phases (main sequence in black, subgiant and red giant 
branch in magenta, central helium burning in blue, post central helium 
burning in red; the pre-main sequence is not shown for clarity). (Bot¬ 
tom) 0.8 M© models computed with Y ini = 0.248, 0.4, and 0.8 (black, 
red, and green, respectively) 


masses are thus expected to give birth to 2P of slightly different 
compositions as a result of the peculiar mixing between individ¬ 
ual winds and pristine material (in other words, a t and a m i n might 
vary for individual polluters). 

Therefore, the abundance variations covered by our grid of 
2P stellar models represent the mean trends for all chemical el¬ 
ements, but do not cover the most extreme abundance variations 
that could be reached in the case of 2P stars formed out of raw 
FRMS ejecta, for example. This can be seen in Fig. [T] where 
the initial chemical composition for our 2P models is compared 
with the maximum abundance variations that have been found 
among main sequence and red giant stars in NGC 6752 (re¬ 
spectively 0.7, 1.7, 1.2, 1.0, 0.46, and 1.8 dex for C, N, O, Na, 


Mg, and Al: ICarretta et al.ll2005i 120071120121) . As discussed in 
iDecressin et al] (l2007bl) . the observed variations in Mg and Al 
require an increase in the 24 Mg burning rate with respect to the 
experimental values (their 60 M© 60rD model), which is not ac¬ 
counted for in the 60 M© 60rC model used in § 13.2.21 However, 
this will not affect any of our findings since the Mg and Al con¬ 
tent has no impact on the stellar evolution tracks or on the stellar 
lifetimes. 

4. Impact of extreme helium enrichment on stellar 
properties all along the evolution: The 0.8M o 
case 

As mentioned above, several studies of the influence of the initial 
helium content on stellar properties can be found in the literature, 
but only for relatively modest helium enrichment or for specific 
evolutionary phases. Here we recall the main trends and explore 
in detail the impact of very high initial helium content (i.e., mass 
fraction > 0.4) all along the evolution. We first describe here 
the 0.8 M© case, which is the typical GC turn-off mass found 
in the literatur^E and later discuss the various evolution paths 
depending on the initial stellar mass (§ 0. The limits of the 
various evolutionary phases are chosen as follows: 

- the ZAMS is defined either by a central mass fraction of hy¬ 
drogen that has decreased by 0.6% compared to its surface 
value, or by a central temperature that has reached 3 x 10 5 6 7 K; 

- the turn-off is reached when the central mass fraction of hy¬ 
drogen is lower than 10 -7 ; 

- the end of the subgiant branch (SGB, base of the RGB) is 
the moment when the mass coordinate at the top of the H- 
burning shell (i.e., when < 10 erg g~ l s~ x ) reaches a local 
minimum. 

- the zero age horizontal branch (ZAHB) corresponds to the 
point where the helium mass fraction at the center of the star 
is 0.8. 

Figure [2] shows the evolution tracks of all the 0.8 M© mod¬ 
els that we computed with the 16 different initial compositions 
described in § 13.2.21 for clarity we also show the tracks of three 
selected cases corresponding to Yi n i=0.248, 0.4, and 0.8. 

4.1. Main sequence 

When the initial helium content is higher, the opacity due to 
Thomson scattering is reduced in the stellar interior, and the 
nuclear burning luminosity increases during the main sequence 
because of the greater mean molecular weighfl. Consequently, 
the evolution tracks are shifted towards higher luminosities and 
effective temperatures, and the main sequence lifetime is signif¬ 
icantly shortened as shown in Fig. [3} while the 0.8 M© model 
with Yi n i=0.248 spends 12.9 Gyr on the main sequence (90 % 
of its total lifetime), those with Yi n i=0.4 and 0.8 have main se¬ 
quence lifetimes of 4.62 Gyr and 158 Myr, respectively, which 
corresponds to 86 % and 60 % of their total evolution. 

During the main sequence, the internal temperature profiles 
steepen and higher central temperatures are reached for higher 
initial helium content (Figs. |4|and0. This modifies the nu¬ 
clear burning mode and the internal stellar configuration, as can 

5 Our 0.8 M© IP model with Y ini = 0.248, Z = 5 x 10 -4 leaves the main 
sequence at an age of 13 Gyr. 

6 The homology relations give L oc p 13 and L oc^F^when the pp chain 

or the CNO cycle dominates; see, e.g jMowlavi et alJ i 1998 ). 
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Fig. 3. Duration of each evolutionary phase of the 0.8 M© models 
as a function of their initial helium mass fraction (main sequence, sub¬ 
giant branch, red giant branch, horizontal branch, and asymptotic giant 
branch are shown in black, magenta, blue, and red, respectively) 



M/M 0 

Fig. 4. Temperature profiles at the zero age main sequence (dashed 
line) and at the turn-off (solid) inside the 0.8 M© models computed with 
the three different initial helium mass fractions 


be seen in the Kippenhahn diagrams presented in Fig. [6] for the 
three 0.8 M© models computed with Yi n i=0.248, 0.4, and 0.8. 
While the first two burn central hydrogen mainly through the 
pp-chain and have a radiative core on the main sequence, the su¬ 
per He-rich model burns hydrogen only through the CNO cycle 
in a convective core (the transition occurs at Y ini = 0.5/0.55). 
In addition, when Yi n i is higher, the effects of central hydrogen 
burning are less pronounced because the initial abundances of all 
the isotopes involved in the CNO-cycle and the NeNa and MgAl 
chains are closer to their equilibrium values. We recall that the 



Fig. 5. Central temperature as a function of the central density along 
the evolution of the 0.8 M© models computed with various initial he¬ 
lium mass fractions from 0.248 to 0.8 (colors correspond to evolution¬ 
ary phases as in Fig.|3; the bend at the end of the RGB is characteristic 
of the helium flash. The straight dotted line corresponds to the limit 
where the pressure given by the non-relativistic perfect gas law is equal 
to the pressure given by the completely degenerate non-relativistic gas, 
and where the mean molecular weight is of the core of a RGB star, i.e., 
log(Tc) = 4.7 + 2/3 log (p c ) 

initial metal mixture at varying Yi n i is obtained from the dilution 
of matter coming from the H-burning regions of massive stars 
with pristine gas. Consequently, the internal chemical profiles 
of the various models strongly differ at the end of the main se¬ 
quence (Fig. [T]), which will have an impact during the advanced 
evolutionary phases. 

4.2. Subgiant and red giant branches 

The models with higher initial helium content cross the 
Hertzsprung gap (subgiant branch) at higher luminosity. Al¬ 
though the first dredge-up occurs in all the 0.8 M© models of 
our grid, its effects on the surface abundances are smoothed out 
for higher Yi n i owing to shallower convective envelopes (see 
in Fig. |7] the vertical arrows that indicate the maximum depth 
reached by the convective envelope during the first dredge-up; 
see also Figs. [6] and [8]) and to initial chemical abundances closer 
to or already at equilibrium. This is quantified in Table [3 which 
summarizes the variations of the surface abundances during this 
phase for the three selected 0.8 M© models. 

The characteristics of the RGB bump are consequently af¬ 
fected. While the 0.8 M© models with Yi n i < 0.4 undergo a drop 
in their total luminosity when the hydrogen-burning shell crosses 
the chemical discontinuity left inside the star by the convective 
envelope during the first dredge-up (Fig. M, the super He-rich 
stars do not and have a constant increase of luminosity along the 
RGB evolution (see the characteristics of the bump in Table [2 
and Fig.©. 

In the same way as for the other evolutionary phases, the life¬ 
times on the subgiant and the red giant branches shorten when 
Yi n i increases (Fig. B- This results from lower degeneracy and 
higher central temperature of the helium core (Fig. ©, which 
lead to smooth central helium ignition occurring at lower stel- 
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Fig. 6. Internal properties (Kippenhahn diagram) of the 0.8 M© models 
with Y ini = 0.248, 0.4, and 0.8 (top, middle, bottom) from the zero age 
main sequence to the end of central-helium burning. Green dashed lines 
delimit hydrogen burning regions (> 10 erg g~ l s -1 ); blue dashed lines 
delimit the helium burning regions (> 10 3 erg g~ l s ~ l ); dotted green 
and blue lines show the maximum energy production of H- and He- 
burning; hatched areas represent the convective zones. Timescales give 
the duration of each evolutionary phase 


lar luminosity for the super He-rich 0.8 M© models, while the 
models with Yi n i below -0.5/0.55 undergo the so-called helium 
flash at the RGB tip (see the different RGB paths in Figsj2]and 
[5]). Quantitatively, it takes 1.18 Gyr (8 % of the total stellar life- 





M/M atar 

Fig. 7. Abundance profiles of the main chemical elements as a function 
of depth in mass at the end of the main sequence for the 0.8 M© models 
with Y ini =0.248, 0.4, and 0.8 (top, middle, bottom). The arrows indicate 
the maximum depth reached by the convective envelope during the first 
dredge-up on the red giant branch 


time) after the turn-off for the 0.8 M© model with Yi n i = 0.248 
to ignite helium, 663 Myr (12 %) when Yi n i = 0.4, and only 
53.8 Myr (20 %) if Yi n i = 0.8. The amount of mass lost through 
the stellar wind during these phases (modeled here according to 
Reimers’ prescription with 77 = 0.5) decreases because of shorter 
timescales and lower luminosity ranging from 0.005 to 0.223 
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M©. Therefore, the total stellar mass at helium ignition is lower 
for lower initial helium content (Fig. □3. 



Fig. 8. Mass layer (in solar mass) within the star reached by convection 
at its maximum extension on the RGB and AGB as a function of the Yi n i 
for the 0.8 M© models 



t (yr) 

Fig. 9. Variations of luminosity during the RGB bump for the 0.8 
M© model with Y ini = 0.248. The luminosity starts dropping when the 
H-burning shell encounters the H-abundance discontinuity left behind 
by the convective envelope during the first dredge-up; it increases again 
when this discontinuity has been crossed. In Table 2 we give the time 
the star spends between the lower and upper horizontal lines, as well as 
the luminosity of the bump indicated by the middle horizontal line 


4.3. Horizontal branch 

The effective temperature and the luminosity of a star on the hor¬ 
izontal branch (HB) directly depend on the total stellar mass, and 
more importantly on the respective masses of the helium core 
and of the envelope above the core at central helium ignition, as 
well as on the opaci t y in the stellar envelope dSalaris & Cassisil 
120051: lMaederll2QQ9l:lKippenhahn et al.ll2013l) . For the reasons 
described above, these quantities strongly vary with the ini¬ 



Fig. 10. Mass of the helium core (black), of the convective envelope 
(red), and total stellar mass (magenta) at the RGB tip (i.e., at helium 
ignition) for the 0.8 M© models as a function of initial helium mass 
fraction. The arrow represents the maximum initial helium mass frac¬ 
tion for the He-flash to occur in the 0.8 M© case 

tial helium content for a given initial stellar mass, as shown in 
Fig.QHIfor the 0.8 M 0 case, with consequences on the HB tracks 
(see Table [2] where we give the effective temperature and corre¬ 
sponding luminosity of the 0.8 M© models on the ZAHB and at 
maximum Teff on the HB). In the same way as for the effective 
temperature on the ZAHB, one can distinguish two regimes. On 
the one hand, the models that undergo the helium flash (Y^i be¬ 
low ~ 0.5/0.55) arrive on the ZAHB with a higher total mass and 
a much lower M core /M env ratio when their initial helium content 
was higher, therefore with a lower T e ff and a higher luminosity. 
On the other hand, for the models that do not undergo the helium 
flash, T e yy(ZAHB) increases with Yi n i because of a higher total 
mass and an increasing core mass at that phase. As a result of 
these two effects, the maximum effective temperature reached on 
the HB by the 0.8 M© models presents a minimum when Yi n i ~ 
0.33. 

In all cases central helium burning occurs in convective con¬ 
ditions (right panels in Fig. li. The resulting CO-core mass in¬ 
creases with increasing Yi n i (Fig. [I2l see also Fig. |TT] for the 
abundance profiles inside the three selected 0.8 M© models at 
the end of central helium burning). 

Finally, the lifetime on the horizontal branch is not very sen¬ 
sitive to the initial helium content in the domain where the stars 
undergo the helium flash (Yi n i below ~ 0.5/0.55), since these 
objects have nearly identical helium-core mass and central tem¬ 
perature at the arrival on the horizontal branch. However, the 
lifetime on the horizontal branch drops for the super helium-rich 
stars. For the 0.8 M© models with Yi n i=0.248, 0.4, and 0.8, this 
phase lasts respectively 94, 91, and 33 Myr (i.e., respectively 
0.7%, 1.7%, and 12.5% of the total stellar lifetime). 

4.4. Asymptotic giant branch 

After central helium exhaustion, all the 0.8 M© models with Yi n i 
below or equal to 0.6 reach the TP-AGB; 0.8 M© models with 


Article number, page 8 of El 

















W. Chantereau et al.: I. Evolution of enhanced helium stars 


He 

(Mass fraction) 

I 

Before 

I 

After 

H 

Before 

[e 

After 

[C/1 

Before 

After 

12 c/' 

Before 

l3 C 

After 

0.248 

0.7514 

0.7396 

0.2481 

0.2581 

0 

0 

90.0 

40.4 

0.4 

0.5994 

0.5905 

0.4001 

0.4090 

-1.30 

-1.36 

5.65 

5.33 

0.8 

0.1994 

0.1994 

0.8001 

0.8001 

-1.61 

-1.61 

3.99 

3.99 


Table 1. Surface abundances of the main elements before and after the first dredge-up on the RGB for the selected 0.8 M© models; the values for 
hydrogen and helium are in mass fraction and those for 12 C/ 13 C are in number fraction. 


Y ini 

log (Lbump/L© ) 

^"bump 

l°g(Ltip/L©) 

log(L/L 0 ) 

ZAHB 

log(T eff ) 
ZAHB (K) 

M to t 

ZAHB (Mq) 

M CO re 

ZAHB (M©) 

log(T eff ) 
max-HB (K) 

0.248 

2.09 

3.9 

3.34 

1.46 

4.19 

0.574 

0.501 

4.21 

0.26 

2.11 

3.6 

3.33 

1.48 

4.17 

0.583 

0.498 

4.18 

0.27 

2.13 

3.4 

3.33 

1.51 

4.16 

0.590 

0.496 

4.17 

0.3 

2.19 

2.6 

3.32 

1.61 

4.10 

0.610 

0.490 

4.10 

0.33 

2.25 

1.9 

3.30 

1.73 

4.02 

0.629 

0.482 

4.07 

0.37 

2.34 

0.99 

3.28 

1.85 

3.90 

0.652 

0.472 

4.09 

0.4 

2.40 

~ 

3.26 

1.92 

3.85 

0.668 

0.463 

4.14 

0.45 

- 

- 

3.21 

2.03 

3.77 

0.693 

0.454 

4.23 

0.5 

- 

- 

3.15 

2.13 

3.75 

0.716 

0.446 

4.34 

0.55 

- 

- 

3.06 

2.24 

3.75 

0.735 

0.443 

4.44 

0.6 

- 

- 

2.95 

2.35 

3.76 

0.745 

0.442 

4.51 

0.65 

- 

- 

2.84 

2.50 

3.88 

0.749 

0.464 

4.57 

0.7 

- 

- 

2.66 

2.65 

4.04 

0.755 

0.490 

4.64 

0.73 

- 

- 

2.60 

2.74 

4.19 

0.758 

0.521 

4.71 

0.77 

- 

- 

2.59 

2.83 

4.41 

0.765 

0.578 

4.69 

0.8 

- 

- 

2.66 

2.76 

4.57 

0.773 

0.627 

4.74 


Table 2. Luminosity and duration (in Myr) of the RGB bump (as defined in Fig. and luminosity at the tip of the RGB. Luminosity, effective 
temperature, total mass, and mass of the core at the ZAHB, and finally maximum effective temperature on the HB for the 0.8 M 0 model as a 
function of Y ini . 


higher Yi n i spend a short amount of time on the early-AGB be¬ 
fore totally losing their envelope. 

During the early-AGB the convective envelope deepens 
again in mass (except in the extreme Yi n i=0.8 case that has ex¬ 
pelled all its envelope at that time; see Fig. [8]). No formal sec¬ 
ond dredge-up occurs, as the convective envelope does not reach 
the regions where helium abundance increases (see the arrows 
in Fig.Qj]), although it barely reaches the region where the CN- 
cycle has operated, and the surface abundances of the carbon and 
nitrogen isotopes slightly change once again. However, this ef¬ 
fect is strongly attenuated when Yi n i increases, as in the case of 
the first dredge-up described above. 

For the 0.8 M© case, the total number of thermal pulses (TPs) 
first increases with the value of the initial helium content owing 
to higher stellar mass at the arrival on the TP-AGB. However 
for Yi n i above 0.45, this is counterbalanced by the greater ratio 
between the mass of the CO core and that of the very shallow 
envelope (Fig. H3, and the number of thermal pulses then de¬ 
creases. Quantitatively, the 0.8 M© with Yi n i=0.248, 0.4, and 
0.6 undergo respectively 8, 18, and 12 thermal pulses, and the 
TP-AGB phase lasts 1.9, 1.35, and 0.07 Myr, respectively. 

Importantly, for Yi n i higher than or equal to 0.65 no ther¬ 
mal pulse occurs and the models spend a very short time on 
the AGB (post-early AGB) and then evolve directly towards the 
white dwarf (WD) stage, which strongly reduces the probability 
of observing super He-rich stars on the AGB phase. 


4.5. White dwarf 

The mass of the degenerate CO-core at the beginning and the 
end of the TP-AGB phase is shown in Fig. [12] for the 0.8 M© 
case. The mass of the CO white dwarf strongly increases with 
increasing initial helium content. The Y = 0.248 model reaches 
the bluest point of the evolution at 14.3 Gyr, and then starts 
the final cooling phase (WD cooling curve) with M core =0.570 
M© whereas for the Y = 0.4 model this happens at 5.41 Gyr 
with M core = 0.659 M© (Fig. fl2l) and at 264 Myr with M cor e = 
0.746 M© for the Y = 0.8 model. Evolution was not computed 
further except for a couple of models, since discussion of the 
WD cooling sequence is beyond the scope of this paper. 

5. Impact of extreme helium enrichment: 
Dependence on the initial stellar mass 

Here we present the combined effects of varying both initial he¬ 
lium and stellar mass over the considered mass domain between 
0.3 and 1.0 M©. We highlight the main points for specific evolu¬ 
tionary phases, and give the corresponding ranges in (M;Yi n i) at 
10 and 13.4 Gyr that delimit the typical ages of Galactic GC£J 
The numbers predicted for (M;Yi n i) at a given age are expected 
to vary with [Fe/H] and are given here for the present grid com¬ 
puted with [Fe/H] =-1.75. The various (M;Yi n i; age) domains 

7 The present models do not include atomic diffusion. Therefore, the 
ages of 10 and 13.4 Gyr should not be directly compared to those de¬ 
rived in the literature for specific GC from models that include the re¬ 
lated processes. Rather, the numbers we give should be compared in 
terms of relative ages. 
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Fig. 11. Abundance profiles as a function of the depth in mass for the 
main chemical elements at the end of the central helium burning for the 
0.8 M© models with Y ini =0.248, 0.4, and 0.8 (from top to bottom). The 
total stellar mass at this evolution point is respectively 0.573, 0.662, and 
0.772 M 0 . The arrows indicate the maximum depth of the convective 
envelope during the early AGB phase 


would also be slightly modified for different prescriptions of the 
mass-loss rates (we note th at in th e presen t computations mass 
loss is treat ed according to lReimerslll975l with Tj = 0.5 on the 
RGB and to lVassiliadis & Wood 1993( ^7the AGB). 

Tracks and lifetime - Similarly to the 0.8 M© case, luminosity 
and effective temperature shifts occur in the Hertzsprung-Russel 



Fig. 12. CO core and total stellar mass at the first thermal pulse, and 
CO core mass after the last one for the 0.8 M 0 models with Y ini between 
0.248 and 0.6 


diagram (HRD) for each given initial mass when Yi n i increases 
(Fig. ns, and the duration of each evolutionary phase is reduced 
accordingly. In particular, the main sequence lifetime strongly 
depends on Yi n i (see the color scale in Fig. [15]). For a canoni¬ 
cal GC age of 13.4 Gyr (see the corresponding isochrone in the 
figure), the stars that are expected to be at the main sequence 
turn-off have (M;Yi n i) between (0.794;0.248) and (0.3;0.719); 
at 10 Gyr, these quantities are between (0.858;0.248) and 
(0.3;0.754). 

RGB bump - The luminosity of the RGB bump increases with 
both initial stellar mass and initial helium content (Fig. ®. 
However, above a certain (M;Y) threshold the stars evolve along 
the RGB without any luminosity drop, and the bump conse¬ 
quently disappears (white area in Fig. [14]) as discussed in § 14.21 
for the 0.8 M© case. The stars that are expected to be at the RGB 
bump have (M;Y ini ) between (0.813; 0.248) and (0.75; 0.295) 
when we consider an age of 13.4 Gyr (see isochrone), and be¬ 
tween (0.881; 0.248) and (0755; 0.336) at 10 Gyr. 

RGB tip and central helium-burning ignition - The (M;Yi n i) 
domain where stars ignite central helium is indicated in Fig. |T5] 
(i.e., in the area found to the right of the curve CO WD / He 
WD curve). For the canonical Yi n i value of 0.248, all the models 
of the considered mass range that ignite central helium-burning 
(i.e., M^i > 0.737 M 0 ) do so in degenerate conditions via the 
so-called helium flash at the tip of the RGB. However, when 
the initial helium increases, low-mass stars evolve faster on the 
RGB; they ignite central helium burning smoothly in less degen¬ 
erate conditions and at a lower luminosity. The limit in terms of 
(M;Yi n i) between the two regimes goes from (0.46; 0.63) to (1.0; 
0.45) for the mass range we are interested in (Fig. HU). 

Position on the horizontal branch - For a given initial stel¬ 
lar mass, the position on the HB strongly depends on Yi n i be¬ 
cause it affects whether or not the star ignites He in degener¬ 
ate conditions, and also depends on the respective masses of 
the helium core and of the stellar envelope when central helium¬ 
burning starts (Sect. 14.31) . In particular, for the models that un¬ 
dergo the He flash, the higher Yi n i, the hotter the position on 


Article number, page 10 ofM 







































































W. Chantereau et al.: I. Evolution of enhanced helium stars 


the HB. At 13.4 Gyr (and at 10 Gyr, respectively), the mod¬ 
els that are in the central helium-burning phases have (M;Yi n i) 
between (0.815;0.248) and (0.670;0.355) (respectively between 
(0.884;0.248) and (0.599;0.458)), and the HB covers a domain 
in T e ff between 9550 and 36310 K (respectively between 5370 
and 36310/37150 K). 

TP-AGB vs failed-AGB - Among the stars that manage to burn 
central helium, a significant percentage finish their life with¬ 
out climbing the AGB. The correspo nding (M;Yj n j) domain of 
these so-called failed-AGB (see, e.g. jGreggio & Renzinil ll990) 
is shown in Fig. Q3] (left of the No AGB / AGB line). These 
stars undergo central helium-burning at very high effective tem¬ 
peratures on the HB; there the stellar envelope surrounding the 
nuclear active core is extremely thin. Consequently, once he¬ 
lium fuel is exhausted at the center, the energy released by the 
contraction of the core and by hydrogen shell-burning mainly 
diffuses across the small, transparent helium-rich external lay¬ 
ers and fails to swell them up. Therefore, the stars first evolve 
at constant effective temperature towards higher luminosity be¬ 
fore moving directly towards the white dwarf region without 
passing through the AGB. This behavior has been proposed as 
a suitable explanati on for the lack of sodiu m-rich stars on the 
AGB in NGC 6752 (ICharbonnel et al.ll2013h . At 13.4 Gyr (and 
at 10 Gyr respectively), the models that are on the AGB have 
(M;Yi n i) between (0.815;0.248) and (0.67;0.355) (respectively 
between (0.884;0.248) and (0.599;0.458)). 

Nature of the white dwarfs - The composition (He vs CO) of 
the stellar remnant strongly depends on the initial (M;Y) as is 
shown in Fig. |T5j In the case of canonical helium abundance 
no He-WDs are expected to be present in GCs today because 
the stellar lifetimes of the corresponding low-mass progenitors 
are much longer than the Hubble tim^E However He-WDs are 
expected to form on much shorter lifetimes for increased initial 
helium abundance above -0.33; in addition, the higher Y^i, the 
lower the initial mass for the star to become a He-WD. The con¬ 
sequences of this behavior on the chemical properties of GC host 
WDs will be investigated in detail in a forthcoming paper. 

6. Conclusions 

This paper presents the first grid of models of low-mass stars 
(between 0.3 and 1 M 0 ) with initial chemical composition de¬ 
rived from the e j ecta of the fast rotating massive star models of 
iDecressin et al.l (l2007albh at [Fe/H]=-1.75 mixed with original 
proto-GC gas. The initial helium mass fraction varies between 
0.248 and 0.8, and the initial abundances of C, N, O, Na, Mg, 
and Al change accordingly. 

The main properties of the models are presented as a func¬ 
tion of initial mass and initial helium abundance. We discuss 
the effects of a very high initial helium content on the evolu¬ 
tion paths in the HRD, on the duration and characteristics of the 
main evolutionary phases, and on the chemical nature of the WD 
remnants. We also give the predicted ranges in initial mass and 
helium content of the stars in different locations of the HRD at 
ages of 10 and 13.4 Gyr that delimit the age range of Galactic 
GCs. 

This grid of models has been successfully used to ex¬ 
plain the lack of Na-rich AGB stars in the GC NGC 6752 

8 For Y ini = 0.248 and Z - 5 x 10 -4 , the maximum initial mass of a 
He-WD progenitor is 0.738 M©, which has a turn-off age of ~ 17.6 Gyr. 
Therefore, for this composition He-WD might be expected only from 
exotic evolution paths; this can happen for example if the low-mass pro¬ 
genitors have undergone extreme mass-loss events , e.g., in close binary 
systems (e.g., Webb inkll 19751 : IStrickler et al.ll2009h . 


(ICharbonnel et al.1120141) . Other aspects related to the properties 
of stellar populations in GCs and of the GC CMDs at various 
phases of their evolution will be investigated in forthcoming pa¬ 
pers in this series. 


Appendix - Content of electronic tables for the grid 

We provide files containing relevant evolution character¬ 
istics (Table [3]) for the 224 models of the present grid. 
All data detailed in Table [3] are available on the website 
http://obswww.unige.ch/Recherche/evol/-Database 
For each mass and initial helium content, models predictions are 
provided from the beginning of the pre-main sequence (along 
the Hayashi track) up to the brightest point of the evolution 
on the AGB. For each model, we have selected 500 points to 
allow a good description of the full raw tracks. First, the most 
important evolution keypoints are determined: 

1. the beginning of the pre-main sequence; 

2. the zero age main-sequence defined as the time when the 
central hydrogen abundance has decreased by 0.003 in mass 
fraction compared to its initial value. 

3. the turning point when the central mass fraction of hydrogen 
is equal to 0.06 on the main sequence; 

4. the end of core H-burning 

5. the bottom of the red giant branch (RGB); 

6. theRGB-tip; 

7. the local minimum of luminosity during central He-burning; 

8. the local maximum of r e ff during central He-burning; 

9. the bottom of AGB: the point with the local minimum of 
luminosity after the loop on HB 

10. the brightest point on the AGB. 

Then, we chose points from the results of evolutionary com¬ 
putation as follows: 

- 100 points are determined alongside the pre-main sequence 
with constant time-steps. The value of each quantity given 
in Table [3] is computed from their values taken from the two 
nearest models calculated with STARE VOL. 

- The same procedure is applied between points 2 and 3 and 
points 3 and 4 where we use step in Xc instead of time; 85 
and 25 points are selected, respectively. 

- Between keypoints 4 and 5 we interpolate 60 points with 
equal time steps. 

- Then 20 points are determined between keypoints 5 and 6 
with equal intervals of log L. 

- 20, 70, and 70 points are selected with equal central hydro¬ 
gen abundance Xq steps between keypoints 6 and 7, key- 
points 7 and 8, and keypoints 8 and 9. 

- 50 points are selected between keypoints 9 and 10 with the 
same log L step on the AGB. 
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Fig. 14. Luminosity of the RGB bump as a function of initial stellar mass and Y ini for those stars that undergo the bump before 14 Gyr (see ages 
given in Gyr along the dotted-line isochrones). The white area corresponds to the (M;Y) domain where stars do not undergo the RGB bump (This 
is also the case for models with an initial mass lower than 0.75 M©) 
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Fig. 15. Nature of the white dwarfs and failed-AGB phenomenon as a function of initial stellar mass and Y ini . The black thin line on the left 
delimits the (M;Y) domains where the stars become either He or CO white dwarfs; it crosses the white line that delimits the domains where helium 
ignition (when it occurs) starts with a flash or in non-degenerate conditions. The curved black thin line separates the stars that climb the AGB from 
those that do not (failed-AGB). Colors depict the age at the main sequence turn-off, with white dashed lines being isochrones with ages indicated 
in Gyr; the white area corresponds to the domain where stellar models predict main sequence lifetimes higher than 13.8 Gyr. 
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Fig. 16. Luminosity of the RGB tip (i.e., when helium ignition occurs) as a function of initial stellar mass and Y ini . The black line delimits 
the domains where helium ignition (when it occurs) starts with a flash or in non-degenerate conditions. The colors indicate the luminosity of the 
RGB tip, with white dashed lines being isochrones with MS turn-off ages indicated in Gyr; the white area corresponds to the domain where stellar 
models predict main sequence lifetimes higher than 13.8 Gyr. 
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Fig. 13. Evolution in the HRD of all the models with initial stellar 
masses between 0.3 and 1.0 M© for Y ini = 0.248, 0.4, and 0.8 (top , 
middle, and bottom, respectively). For the very low-mass models, we 
stop the evolution if it exceeds 14 Gyr. At this age, the (0.3, 0.248) and 
(0.3; 0.4) models are still on the pre-main sequence, therefore they are 
not present in the HRD. 
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